Critical behaviour of driven bilayer systems: 
A field-theoretic renormalisation group study 
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Abstract. We investigate the static and dynamic critical behaviour of a uniformly 
driven bilayer Ising lattice gas at half filling. Depending on the strength of the 
interlayer coupling J, phase separation occurs across or within the two layers. The 
former transitions are controlled by the universality class of model A (corresponding 
to an Ising model with Glauber dynamics), with upper critical dimension dc — 4. The 
latter transitions are dominated by the universality class of the standard (single-layer) 
driven Ising lattice gas, with dc = 5 and a non-classical anisotropy exponent. These 
two distinct critical lines meet at a non-equilibrium bicritical point which also falls 
into the driven Ising class. At all transitions, novel couplings and dangerous irrelevant 
operators determine corrections to scaling. 
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Introduction. Driven diffusive lattice gases (DDLG), introduced by Katz et al. 
||T|] to investigate far- from- equilibrium properties of interacting many-particle systems, 
are deceptively simple generalisations of the familiar equilibrium Ising lattice gas 
Particles diffuse on a lattice, controlled by not only the usual interparticle attraction 
and thermal bath (at temperature T), but also a uniform external force. In conjunction 
with periodic boundary conditions, the latter drives the system into a non-equilibrium 
steady state with non-trivial current. A remarkable range of novel collective phenomena 
emerges, some aspects of which are now well understood while many others remain 
mysterious 0]. In particular, when further 'slight' modifications or generalisations are 
introduced, simulations often contradict equilibrium-based expectations. In this sense, 
an intuitive understanding of non-equilibrium steady states is still lacking. 

While most simulations of the original model were carried out on a single two- 
dimensional {d = 2) square lattice, Monte Carlo studies have been reported recently for 
two coupled driven Ising lattices stacked to form a bilayer structure 0, ||, Q . Upon tuning 
the (Ising) inter-layer coupling, the system makes a transition from a homogeneous 
high temperature phase into two distinct ordered states at low temperatures. A simple 
phase diagram was found [^, displaying two lines of continuous transitions which meet 
a line of 'first-order' transitions at a 'bicritical' point. The critical properties associated 
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with the two hnes were conjectured. In the most recent study [^], anisotropic intra- 
layer couphngs were introduced, and certain critical properties were measured. In this 
letter, we report results from a field-theoretic renormalisation group (RG) study for the 
continuous transitions of this model. Even though we essentially confirm the original 
conjecture ||^, some novel and curious features emerge near the bicritical point and along 
the line associated with repulsive inter-layer couplings. After a brief description of both 
the microscopic model and the field-theoretic description, we present our results. 

In all simulation studies, the 'microscopic' model consists of two fully periodic 
Li X L2 square lattices, arranged in a bilayer structure (effectively, an Li x L2 x 2 
system). The sites, labeled by (ji,j2,j3), with ji^2 = l,--,-Z^i,2 and js = 1,2, may be 
empty or occupied. Thus, the set of occupation numbers {n{ji, j2, j3)}, where n = 
or 1, completely specifies a configuration. To access critical points, we use half-filled 
systems, i.e., = -^1-^2- The particles interact, so that the Hamiltonian is given by 
Ti, = —jQ^nn' — J^nn" , where n and n' are nearest neighbours within a given layer, 
while n and n" differ only in the layer index. All studies focus on attractive intra-layer 
interactions, Jq > 0, while J can be of either sign. Since intra-layer anisotropies generate 
no qualitatively new features, Jo may be set to unity. The equilibrium phase diagram 
in the J-T plane is easily obtained, with some exactly known features. For example, a 
second-order transition occurs at Tq — 0.5673/A;b and J = |]^. 

To access the most interesting non- equilibrium properties, a conserved dynamics 
must be imposed. Typically, Kawasaki spin exchange with Metropolis rates is employed, 
i.e., particles hop to nearest neighbour holes with probability min{l,exp {—AH/kBT)}, 
where A7i is the energy change associated with the move. To model the effects of the 
drive, we add ±Eo to AH for hops against /along, say, the 1-axis [Q, interpreting the 
particles as 'charged' in the presence of an external 'electric' field {Eo, 0,0). Note that, 
with Metropolis rates, it is possible to study the 'infinite' case: jumps against the 
field are simply never executed. When so driven, the phase diagram in the T-J plane can 
be found (schematically shown in Fig. 1 ^, ||). At high T, the system is disordered 
(D). At low T, the system phase separates: For sufficiently repulsive J, the fully ordered 
(T = 0) state displays densities 1 and in the two layers, so that this phase is labeled 
'full-empty' (FE). On the other hand, for J > and low T, each layer phase separates 
individually, resulting in strips of particles (of width L2/2, at T = 0) 'on top of each 
other', and aligned with the drive. Thus, we label this state the 'strip' (S) phase. In 
the following, we invoke field-theoretic renormalisation group techniques, to investigate 
universal properties associated with the critical points. 

Model equations. This analysis has already been initiated in p, ^. We define 
the single-layer magnetisations ifi{x), i = 1,2 as the coarse-grained versions of 
2n{ji, j2,i) — 1, with the in-layer coordinate x generalized to d dimensions. To insure 
the proper Eq = limit, we construct a Landau-Ginzburg- Wilson (LGW) Hamiltonian, 
Tic, which contains all terms, up to fourth order in cpi and second order in Vy^i, 
compatible with stability requirements and symmetries of the microscopic model. 
As usual, the explicit relationships |§] between the coarse-grained couplings and the 
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Figure 1. Schematic phase diagram in the J,T plane, for Eo = oo. The D-FE (D- 
S) transition hne is shown solid (dotted). The S and FE phases are separated by a 
(dashed) first order line. The junction of the three lines marks the bicritical point. 



microscopic parameters (J, Jo,T) are not needed. The next step is to incorporate the 
dynamics and the (coarse-grained) drive E (oc tanhE'o), following Ref. [Q. The result 
is a set of Langevin equations for the fields (pi{x,t). Focusing on the deterministic 
evolution first, the basic equation for (pi (with a similar one for ip2) reads dtfi = 
f (1 - ip2) [5Hcl5ip2 - 5nc/5ipi]+-i [V^ {SHc/Sipi) + V2Edip]\ . Here, the first term 
models the energetics of inter-plane jumps, with a relaxation constant denoted by A. 
The second term reflects the (conserved) in-plane dynamics, being the coarse-grained 
version of the usual (single-layer) DDLG, with diffusion coefficient 7. The symbol d 
indicates a spatial gradient along the direction of the drive. Some numerical factors 
appear for later convenience. Anticipating the very different roles played by the total 
and 'staggered' magnetisations, we introduce S(x, t) = [ipi{x,t) + ip2{.x,t)\/ \/2 and 
A = — Lp2\l^2. Finally, we incoporate the appropriate Langevin noise terms into 
the equations of motion, so that the starting point of our field theoretic analysis is: 




= -A (rA - V" ) A + ^ A^ + ^ A 
+ A(EA9S + E'S9Aj +Ca 
with stochastic white noise characterised by (Ce) = = (Ca), and correlations 

(Ce(x, t) Ce(x , t')) = -27 V b{x - f ) b{t - t') , 
(Ca(x, t) Ca(x , t')) = 2A 6{x - x) 6{t - t') . 



(1) 



(2) 



(3) 
(4) 



Some comments are in order. First, E = E = E' = E in mean-field, but suffer different 
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renormalisations, so that our parameter space must be generalized. Second, for zero 
drive, the D-FE (D-S) transitions are marked by the vanishing of (^^e); i-e., A (S) 
becoming critical. Third, the local E is conserved and unaffected by inter- layer jumps; 
hence, its relaxation is diffusive and independent of A. Indeed, for u = and vanishing 
drive, Eq. (|I|) reduces to that in model B dynamics |^. Now, for non-zero drive, the 
terms containing only S-fields are precisely the standard Langevin equation for the 
DDLG |lO|, |ll|]. Therefore, Eq. is supplemented with a conserved stochastic noise 



term. On the other hand, for ^ = and zero drive (0) reduces to the purely dissipative 
dynamics of model A [^, controlled only by the inter-plane jump rate A. Thus, non- 
conserved white noise has been added. The cross correlations (Cs Ca) vanish, as a 
consequence of the statistical independence of intra- and inter-layer jumps. Finally, the 
drive will induce the standard DDLG spatial anisotropics, to be discussed next. 

The D-S critical line. We begin with the D-S transitions, where the conserved 
field E becomes massless, while A remains non-critical. In the presence of the driving 
term, the characteristic anisotropics of DDLG's emerge. The system softens only in the 
spatial sector transverse to the drive, so that transverse and longitudinal momenta scale 
differently: q\\ ~ q]^'^, where k is known as the anisotropy exponent (k = 1 in mean- 
field) [[10| , pJ| . Consequently many longitudinal non-linear terms become irrelevant. 
Next, since the non-conserved field A is massive, loop diagrams containing internal A 
propagators are m/rared- convergent and may be ignored. Thus g and, remarkably, (a 
may be neglected. Alternatively, the same conclusions can be reached by focusing on 
ultraviolet singularities in perturbation theory (see Ref. [0). The result is that, within 



the E sector, the analysis is identical to the standard single-layer DDLG case [TO, [TI 
These considerations can be summarized in the effective theory, 

9tE = 7 [c 9^ + Vi (r^ - Vl)] E + 7 I ^E^ + , (5) 
atA = -A(rA-V') A-a|ae2 + AEA9E + Ca , (6) 



where c > is introduced to account for anomalous anisotropy jT^, |rT|. The D-S 
line itself is associated with a vanishing (renormalized) r±, while ta > (so V^A is 
also irrelevant). Denoting a momentum scale by fi (~ q±), a consistent set of naive 
dimensions emerges: 7 ~ u ~ fi^, A ~ /i^, and c ~ /i". As in the standard 
case, the most relevant coupling is £^ = E'^jc^l'^ (~ yU^~'^ naively), so that the upper 
critical dimension is = 5. In addition to c, only the coupling o = go/ EE exhibits a 
non-trivial flow under the RG. The remaining effective non-linearities all have lower 
scaling dimension: u,u, EE, E'E ~ and EE' ~ Of course, the static 

coupling u drives the phase transition and therefore constitutes a dangerously irrelevant 
variable p^. 

Remarkably, the analysis can be performed to all orders in e = — d. In particular, 
the E-loops which modify the new vertex (E) are the same as those affecting c. Thus, the 
same power series, p{S), enters the RG flow equations for both renormalised couplings: 
(3£ = ixd^E = —£[e + \p{£)] and (3^ = /xS^a = a(l + a)p{£). For ci < 5, the stable RG 
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fixed point turns out to be p{S*) = —2e/3 and cr* = —1 to all orders, whereas for d > 5, 
both couplings tend to the Gaussian fixed point 0. Notice, however, that the novel 
bilayer coupling a does not enter any singular diagram for the two-point functions. 
Consequently, despite its non-trivial fixed-point value, it does not affect the scaling 
exponents along the critical E line, which are just those of the standard DDLG. Even 
for d < 5, the only non-trivial exponent is n. The two-point correlation function for 
the S fields scales as Css(g', a;, r_L) = qj^ Cj^ Qw/Qx''^, "^/^l); corresponding 

to the transverse critical exponents ri± = 0, = 1/2, and z± =4. As a consequence 



of Galilean invariance |jTO|, the exponent k is fixed by a scaling relation to 1 + e/3 for 
2 < < 5 (k = 1 for > 5). 

The bicritical point. Remarkably, these features remain valid at the point where 
both fields are critical: r/^ = r± = 0. Formally, this follows from the observation that the 
above scaling dimensions still hold, and no additional diagrams appear. The conserved 
field S essentially slaves the non-conserved field A. The larger critical dimension d^ = 5 
still dominates the scaling behaviour, keeping the static couplings u, u, g, and g irrelevant 
(with g appearing merely through a). Thus, there are no non-trivial renormalisations in 
the two-point function for the A-fields, so that its scaling is just isotropic and mean-field- 
like, C/^A{q,uj,rj^) = (^qr^^^^'^ , uj/q^^. This corresponds to a non-conserved 
Gaussian theory: 77 = 0, z/ = 1/2, and z = 2 near five dimensions. 

The D-FE critical line. Finally, we turn to the case where only the A field is 
critical: — > 0, > 0. Though E remains massive, it is a conserved field and hence 
a slow variable. This situation is reminiscent of models E or G in equilibrium critical 
dynamics [^j. The key difference resides, however, in the external drive rendering the 
system far-from-equilibrium and manifestly anisotropic. At the same time, the coupling 
to the diffusive mode sets us outside the framework of the simple nonequilibrium kinetic 



Ising models investigated in |T^. So, a full RG calculation is required to determine 
whether model-A relaxational kinetics still persists. 

With E non-critical, we may safely neglect the terms V^E and V^E^. Naively, 
scaling is isotropic along this critical line, so that {q\\,q±) ^ ji, uj ^ /i^, and 
(A, 7) ~ The dominant non-linear coupling is now g ~ with an upper critical 
dimension d^ = 4. Power counting yields u ~ /i^"*^ which therefore becomes irrelevant. 
However, beyond g, several additional marginal effective couplings appear, consisting 
of combinations such as EE, EE, E'E, gE/E. Hence the effective critical theory still 
contains two anisotropic propagators and six non-linear vertices: 

9,E = 7(c92 + V1)e + 79(^|e^ + |a2)+Ce, (7) 
dtA = -X(rA-ad'-Vl)A-x(^A^ + ^A T? 



.6 2 

+ A(EA(9E + E'EaA) +Ca , (8) 

Similar to the D-S line, all diagrams that contain the conserved noise Cs are non- 
singular. Thus, in contrast to, e.g., model C, becomes irrelevant here. Instead, the 



Critical behaviour of driven bilayer systems 



6 



dynamics of the conserved field is dominated, via the couphng E, by the fluctuating, 
non-conserved, critical A. This remarkable feature may be observable in simulations of 
the E-S correlations, which should develop anomalous (beyond the generic variety 
long-range components. 

Carrying out an e (= 4 — (i)-expansion at the one-loop level, we find 35 non-trivial 
Feynman graphs, complicated further by the two distinct anisotropies in A and S. To 
be brief, we only report the salient features and leave details to be published elsewhere 



15| . The RG fixed points are determined from seven coupled RG flow equations for five 



marginal non-linearities {g = g/a^^"^, f = EE/a^^"^, f = EE/a^^"^, f = E'E/a^/'^ and 
h = gE / Ea^^"^) and two dimensionless ratios (v = c/a, w = 'j/ \). All equations involve 
(f , w) in algebraic expressions and, in general, fixed points (f *, w*) are not simply O (e). 
Aside from the Gaussian fixed point, the most interesting non-trivial zeroes of the RG 
/^-functions are the symmetric fixed points (/* = /'* = (fh/gv)* ^ 0) and the fixed 
points restoring Galilean invariance {h* = 0, /'* = w*f*). Unfortunately, none are 
stable under the RG flow. The completely stable fixed points we found are essentially 
in the class of model A: g* = e/3 + O(e^) and /* = /'* = h* = 0. Remarkably, two 
unusual features emerge: First, {v*,w*) remain undetermined by the flow equations. 
Instead, they are constrained by the stability analysis alone, resulting in a stable fixed 
domain V Second, for each {v*,w*) G V, the coupling /* assumes a non-trivial 

-|- O(e^). However, none of these novel aspects 



3i«*2 3w* 



value e [(^1+^. ^.+^t+^)2 (i+v^)2_ 
affect the leading exponents, as we summarize below. 

Aside from technicalities, the features of the D-FE and D-S lines are similar: The 
critical exponents are determined by the critical field alone, while the noise of the massive 
mode becomes irrelevant. A single coupling survives, through which the critical field 
influences the dynamics of the massive mode, but not vice versa. Thus, on the D-FE 
line, we find standard model A behaviour: Caa(i?, t^, = q^^^'^^^ Ca {q'^^ ■, ^l^f)-, 
with ?7 = 0, v'^ = 2 — e/3 and 2; = 2, up to 0(e) (for 2 < c? < 4) and trivial 
anisotropic properties. Meanwhile, the leading singularities of the S field are of the 
simple diffusive type: C-£Y;{q,u!) = gj^ C^^u/q"^) (with no anomalous anisotropy: 
K = 0). We stress, however, that unlike the usual kinetic Ising model, our bi-layer 



system exhibits non-vanishing and anomalous three-point correlations |]T6| — an obvious 



signature of the external drive. These should be clearly detectable in simulations. In 
addition, corrections to scaling are expected to be distinct from those of model A. 

Summary and discussion. To conclude, we find that the universal critical behaviour 
of a driven, bilayer Ising lattice gas is controlled by two distinct universality classes, 
depending on whether the system orders into S or FE configurations, i.e., whether 
the conserved sum, E, or the non-conserved difference. A, of the coarse-grained layer 
magnetisations becomes critical. In both cases, the fixed point theory shows the critical 
field decoupling completely from the non-critical one, as the noise of the latter becomes 
irrelevant. Thus, all along the D-S line, the critical dimension is d^ = 5, and the 
scaling exponents for the critical E-field are just those of the "standard" DDLG. In 
contrast, for the D-FE transitions, the leading singularities are that of model A for the 
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A-field, with = 4 and isotropic exponents. At the bicritical point, though both S 
and A are critical, the larger critical dimension is df = 5, so that the system again 
displays DDLG-like critical properties. However, at all critical points, a single non- 
trivial operator survives which provides a one-way coupling of the critical field into the 
dynamics of the non-critical one. For the D-S transitions, this operator is a = gc/EE, 
while its counterpart at the D-FE transitions is / = EE/a^/'^. Interestingly, both have 
negative fixed point values, indicating that the three-point couplings do not all share the 
same sign. The stability domains of these fixed points are bounded, so that trajectories 
starting outside the stable region can run off to infinity |jl5|. Even for stable trajectories, 
the presence of numerous other, unstable fixed points tends to enlarge crossover regimes, 
shrouding the asymptotic universal scaling behaviour. This may in fact be at the core 
of the inconsistent exponent values observed in Ref. A better understanding of 
the bare couplings, associated with the microscopic Hamiltonian, would be extremely 
helpful. 
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